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Context. Sun-like stars are promising candidates to host exoplanets and are often included in exoplanet surveys by radial velocity 

(RV) and direct imaging. In this paper we report on the detection of a stellar companion to the nearby solar- analogue star HD 104304, 

which previously was considered to host a planetary mass or brown dwarf companion. 

Aims. We searched for close stellar and substellar companions around extrasolar planet host stars with high angular resolution imaging 

to characterize planet formation environments. 

/Wefhods. The detection of the stellar companion was achieved by high angular resolution measurements, using the "Lucky Imaging" 

technique at the ESO NTT 3.5m with the AstraLux Sur instrument. We combined the results with VLT/NACO archive data, where 

the companion could also be detected. The results were compared to precise RV measurements of HD 104304, obtained at the Lick 

and Keck observatories from 2001-2010. 

Results. We confirmed common proper motion of the binary system. A spectral type of M4V of the companion and a mass of 0.2 1 Mq 

was derived. Due to comparison of the data with RV measurements of the unconfirmed planet candidate listed in the Extrasolar Planets 

Encyclopaedia, we suggest that the discovered companion is the origin of the RV trend and that the inclination of the orbit of / a 35° 

explains the relatively small RV signal. 
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1. Introduction 

HD 104304 (LHS 5206, HIP 58576) i s a high proper motion 
star at a distance of 12.91 ± O.lSpc (Ferryman et al. 1997|l. 



Because of its proximity, brightness, and observability from both 
hemispheres, it has been the subject of numerous studies. It is a 
chromosph erically inactive star with log R'(HK) = -4.933 (Gray 
|et al.|2006| and with an X-ray luminosity of Lx = 1.2 X lO^'^ 
( [Schmitt & Liefke|2004| l. It has a spectral type of G8IV to KOI V, 
log g = 4.35 to 4.15, and [Fe/H] = 0.22 to 0. 18 (IGray et al.|2006' 



CenaiTO et al. 2007), and is a thin disk member (" Soubiran &| 
|Girard.2005). With an age of ^8.5 Gyr and a mass of =sl.01 Mq 



( Takeda et al.|2007 1 it is considered a close analogue to the Sun 
at approximately twice its present age, and is a target of a high- 
precision Doppler survey aimed at identifying planetary com- 
panions to Sun-like stars ( [Wright et aL]|2()D4) l. HD 104304 has 
also been suggested as a target to search for habitable planets 
with future space missions ( [Lopez et al.|2005| l. 

In 2008 we started a high-angular resolution imaging survey 
of exoplanet hosts stars with the goal to identify close stellar 



and substellar companions. The first results of this survey were 
reported by|Daem gen et al.[ ( |2009| . HD 104304 was included in 
the survey because of the identification of a potential substellar 
companion by Ma rcy|p007 1. 

In the following we will use the designation HD 104304A 
for the primary star and HD 104304B for the companion. 



2. Observations and data reduction 

2.1. AstraLux Sur 

Observations at the ESO 3.5m NTT were obtained with the 
'Lucky Imaging' instrument AstraLux Sur, a high-speed electron 
multiplying camera for Lucky Imaging observations (Hippler 
|etal.|2009 ). The instrument is a clone of the common user instru- 
ment AstraLux Norte at the Calar Alto 2.2 m telescope ([Hormuth 
et al. 2008 ). It provides a field of view of 16 x 16 arcseconds. 



* Based on observations made with ESO Telescopes at the La Silla 
and Paranal Observatory under programme IDs 083.C-0145 and 084.C- 
0812, and on data obtained from the ESO Science Archive Facility. 



HD 104304 was observed as part of a larger sample of (can- 
didate) exoplanet host stars stars on April 1 1, 2009, and again on 
February 6, 2010 in i' and z'-band. Total integration times were 
30 s corresponding to the best 10 % of 10000 frames with an in- 
dividual texp - 0.030s. Figure 1 shows the AstraLux images of 
the binary as well as the VLT/NACO data (see section 2.2). 
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Fig. 1. AstraLux z' and NACO/SDI images of the binary using a logarithmic intensity scale. The field of view is 3" x 3", and North 
is up and East is to the left. The linear projected separation of the components is 12.8 AU, 14.2 AU and 14.8 AU from the left to the 
right image. 



The data were reduced by the AstraLux online data reduc- 
tion pipeline as explained by |Hormuth et al.| ( |2008| l. Pixel scale 
and image orientation were determined by observations of star 
clusters with well-determined astrometry. 

2.2. VLT/NACO 

After the detection of the companion, we searched the VLT 
archive for high-angular resolution observations of HD 104304. 
VLT observations of HD 104304 using the Spectral Differential 
Imager (SDI) mode of NACO were obtained on January 18, 
2006 in service mode under the programme-ID 076.C-0762(A) 
for Damien Segransan, and retrieved from the ESO archive. 
NACO/SDQ provides four simultaneous images in three differ- 
ent filters in the near infrared H-band. Total integration times 
were 20x172.5 s. The companion was clearly detected in the 
NACO images. 

We also checked the deeper NACO/SDI frames for addi- 
tional closer and fainter companions by reducing the data with a 
dedicated SDI+ADI pipeline in IDL, see Figurel2] For details on 
this pipeline see Janson et al. (2007). These data would enable 



us to detect a companion 7.5 mag fainter in H than HD 104304A 
with a detection limit of 3 sigma at an angular separation of 0.3" 
(as assumed for the unconfirmed planet in http://exoplanet.eu). 
No evidence for additional companions could be found. 

2.3. Determination of binary parameters 

FigurefTlshows a comparison of the AstraLux z', and the NACO 
images of the HD 104304 binary. In the NACO data set, we 
analysed the four acquisition images with the shortest inte- 
gration time to derive precise relative astrometry and photom- 
etry for both components. For the analysis a pixel scale of 
17.32 mas/pixel was assumed. Aperture photometry was con- 
ducted to determine the brightness ratio. The position of com- 
panion and primary star were determined by fitting a Gaussian 
profile. 

In the AstraLux data set, the procedure to determine the 
brightness ratio was slightly different, as the companion is situ- 
ated inside the seeing halo of the parent star Aperture photome- 
try was conducted, taking also into account the halo background. 
This background was determined on the opposite site of the pri- 




For more information on NACO/SDI 

[http://www.eso.org/sci/facili ties/paranal/instruments/naco/inst/| 



Fig. 2. With the SDIh-ADI pipeline reduced NACO/SDI frames 
of HD 104304. In the center of the image the residuals of the 
subtracted primary star can be seen. The white and black dots on 
the right half side of the image are the companion on the different 
orientated images. 



mary star from where the companion was situated, assuming a 
symmetric Point Spread Function (PSF). To determine position 
angle and separation, PSF subtraction was applied to the primary 
star The PSF for subtraction was created by mirroring the half 
of the primary star without companion to the opposite side. On 
this PSF subtracted image, the position of the binary companion 
was determined by fitting a Gaussian profile. The position of the 
parent star was determined by fitting a Gaussian profile to the 
original data. The results can be found in Table [T] 

3. Physical properties of tlie HD 104304 binary 

3.1. Common proper motion 

HD 104304 is a high proper-motion star with ;Ura - 141.98 + 
0.71 mas/yr and //dec = -483.34 ± 0.42 mas/yr. In the 3.25 years 
that passed between the VLT/NACO and the first AstraLux ob- 
servations, HD 104304 moved 461.44 + 2.3 mas to the East, and 
1570.9 ± 1 .4 mas to the South. In the same period, the separation 
between HD 104304A and B changed from 995+2 mas to 1 104± 
2 mas, and the position angle increased from 300.19 ± 0.01 deg 
to 333.02 + 0.08 deg (see TablefT). If the companion had been a 
background object, the corresponding numbers would have been 
a change in separation to 1404 + 2 mas and in position angle 
to 314.47 + 0.03 deg, which can be exluded with high signifi- 
cance. Thus both sources form a physical binary, and the change 
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Table 1. Binary properties 



Date 


Filter 


separation 


PA 


mag. diff. 






["] 


[deg] 


[mag] 


2006-01-18 


NACO SDI 


0.993 ± 0.005 


300.1 ±0.1 


3.80 ± 0.05 


2009-04-11' 


SDSS i' 
SDSS z' 


1.104 ±0.002 


333.02 ± 0.08 


5.5 ±0.3 
4.8 ±0.2 


20 10-02-06 '-2 


SDSS z' 


1.15 ±0.02 


340.45 ±0.15 





' Separation and position angle are the mean values of these parameters in i' and z'-band. 

^ The magnitude difference of this data set was not determined as the seeing was much worse than in the 2009 data. 



in position angle can be attributed to orbital motion around the 
common center of mass of both components. 

3.2. Spectral types and component masses 



As stated by |Kraus & Hillenbrand ( 2007[ l, the literature on in- 
trinsic magnitudes and basic astrophysical parameters of main- 
sequence stars comprises a rather heterogeneous set. In order to 
minimize systematic errors in the determination of the properties 
of the HD 104304 binary system, we base the following analy- 
sis therefore on the homogenized set of parameters as compiled 
and computed by Kraus & Hillenbrand] ( 2007 1. Assuming that 
the average of the brightness differences in the three NACO/SDI 
spectral bands is close to the brightness difference in H-band, 
we can derive the absolute H-band magnitudes of both compo- 
nents based on the 2MASS photometry of the unresolved binary 
star. The absolute H-band magnitude of HD 104304B derived 
this way is consistent with a spectral type of M4V. 

To determine the i' magnitude of the binary companion the 
magnitudes given by |Just & JahreiB] (J2008jl of the unresolved 
HD 104304 binary were used. By a comparison of the magni- 
tude with spectral types given by Kraus & Hillenbrand' (2007") 
the spectral type of the companion was derived. Component A 
has a spectral type of G8IV to KOIV which leads also to a spec- 
tral type of M4V of HD 104304B in i'-band. We note that for 
mid M-dwarfs, i'-band magnitudes vary strongly with spectral 
type. Thus an uncertainty by -1-0.4 mag in i' corresponds to an 
uncertainty by less than half a spectral-subclass. By interpolat- 
ing the mass scale in Table 5 from Kraus & Hillenbrand ( 2007| l 
we derive a mass of 0.21+°J» Mq for HD 104304B. 

The derived absolute magnitudes of both components in H 
and i'-band are listed in Table |2] 

Table 2. Component Photometry (absolute magnitudes) 



Component 



niH 
[mag] 



[mag] 



SpT 



A 
B 



3.5: 
7.3: 



0.2 
0.2 



4.6 d 

10.1: 



0.3 

:0.4 



G8IV to KOIV 

M4V 



1.01 Mo 

0.21!°ii^Mo 



3.3. Orbital parameters 

While the projected separation of the binary changed only by a 
relatively small amount from 995 mas to 1 104 mas, the position 
angle increased quite substantially by 32.83 deg over a period of 
3.25 years. This gives some evidence that the orientation of the 
orbit is closer to face-on than to edge-on (i.e. / <45°). 

The HIPPARCOS p arallax of HD 104304 is 77.48 ±0.80 mas 
( [Perryman et al.|1997| l, corresponding to a distance of 12.91 ± 



0.13 pc. Hence the angular separations correspond to projected 
separations of ^12.8 AU to s!l4.8 AU. 

Taking into account only the change of the position angle 
during the 3.25 years, assuming a circular orbit and a semimajor 
axis of a;14AU, we would expect an orbital period of a;36yr. 
With these assumed orbital parameters, Kepler's third law yields 
a system mass of ^2.1 Mq. Hence, the real eccentric orbit has 
a smaller semimajor axis or the orbital period is larger, as the 
system has a mass of x!l.2M0 (see Tablel2]). 

To precisely determine the orbital parameters, more astro- 
metric data points are necessary. Therefore astrometric monitor- 
ing of the binary over the next years is required. 

3.4. Comparison with data of unconfirmed planet candidate 

Based on its RV trend, HD 104304 has been identified as a can- 
didate host of a substellar object (Marcy 2007 ). This RV trend 
was measured during 5 years, from 2001-2006. In these 5 years 
a change of a;500 m/s in the RV signal was detected. 

We calculated the estimated radial velocity curves for a bi- 
nary system with 1 .22 Mq, different orbital periods between P = 
35 years and P = 60 years different inclinations from / = 5 ° to 
/ = 50 ° and different eccentricities . 

To these theoretical RV curves, we compared precise RV 
measurements of HD 104304, taken over a period of 9 years. RV 
measurements at the Lick observatory have been obtained dur- 
ing a period of 6.5 years, from May 2001 to December 2007. At 
the Keck observatory, RV measurements have been taken over 
a period of 3 years, from December 2006 to January 2010. As 
can be seen on Figure l3] the calculated RV curve for an inclina- 
tion of / = 35 °, an eccentricity of e = 0.29 and a period of P = 
48.5 years compares well with the RV data. The identification of 
the unconfirmed planet by Marcy ("2007) was also based on these 
RV measurements from the Lick observatory, only on a shorter 
timescale than the RV data presented here. Therefore we sug- 
"gest that the extrasolar_planet candidate listed in the Extrasolar 
Planets Encyclopaedicj^ is identical with the stellar companion 
"presented in this paper. 

While a plausible set of orbital parameters is derived, the 
analysis also hints that thus far only around 20% of the orbit has 
been covered by RV monitoring, and less than 10% of the orbit 
by astrometric observations. Thus further monitoring is required 
in order to derive a robust set of orbital parameters. 

Theoretical calculations were made in order to check if a real 
planet/brown dwarf companion orbiting HD 104304 would have 
been visible in the NACO/SDI frames. 

The reduced NACO/SDI data (see Figurel2]i would enable us 
to detect a companion 7.5 mag fainter in H than HD 104304A 



^ The Extrasolar Planets Encyclopaedia can be found at 
http://exoplanet.eu and is maintained by Jean Schneider. 
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Fig. 3. Left: Theoretical RV curve of the binary system with a period of P = 48.5 years, an inclination of / = 35° and an eccentricity 
of e = 0.29. Triangles are the Lick RV measurements from 2001-2007 and asterisks are the Keck RV measurements from 2006- 
2010. Right: Residuals of the model. Diamonds are the residuals of the model with respect to the Lick measurements, asterisks with 
respect to the Keck measurements. 



at an angular separation of 0.3", i.e. a source with an absolute 
magnitude mn = 11.0 mag for a detection limit of 3 sigma. In 
the age range 5 to 10 Gyr, according to Baraffe et al. (2003), this 
corresponds to a star with 83 Mj^p (0.08 Mq). Thus while we can 
exclude an additional stellar companion with a mass >83 Mjup 
at a projected separation >3.9 AU, a substellar (brown dwarf or 
planetary mass companion) would not be visible in the NACO 
images. 



4. Conclusions 

The combination of NACO/SDI archival data and AstraLux Sur 
data made it possible to confirm the stellar companion around 
HD 104304. Based on the photometry, a spectral type of M4V 
and a mass of 0.21 Mq were derived for HD 104304B. 

The primary star is well studied, e.g. it has a known metal- 
licity of [Fe/H] = 0.22 to 0.18 (|Gray et al. '20061 [Ceharro et al. 
2007| l, an age of ^ 8.5 Gyr ( [TakedaeTa l. 2007 ) and a mass nearly 
equal to that of our Sun. This knowledge is a valuable asset for 
the further characterisation of any stellar or substellar compan- 
ion in this system. 

By modeling the theoretical RV curve and comparing it to 
the available RV observations obtained over a period of 9 years, 
we have shown that the observed RV trend could be explained by 
the detected stellar companion. A theoretical RV curve with an 
eccentricity of e = 0.29, an orbital period of P = 48.5 years and 
an inclination of the orbit of / = 35° fits best the RV measure- 
ments. This gives evidence that the planet candidate listed in the 
Extrasolar Planets Encyclopaedia around HD 104304 is identical 
with the presented stellar companion. Astrometric monitoring of 
the binary system over the next few years is necessary to con- 
firm this statement and to precisely determine the orbit and its 
orientation. 

In the case of confirmation of a planet around HD 104304, 
this system would be highly interesting in the context of planet 
formation and migration around close binaries with solar ana- 
logue host stars. 



AstraLux Sur as well as all the technical staff at ESO La Silla and MPIA involved 
in the project. 
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